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ABSTRACT 

We present an analytical approximation formula for the growth function in a spatially 
flat cosmology with dust and a cosmological constant. Our approximate formula is 
written simply in terms of a rational function. We also show the approximate formula 
in a dust cosmology without a cosmological constant, directly as a function of the 
scale factor in terms of a rational function. The single rational function applies for 
all, open, closed and flat universes. Our results involve no elliptic functions, and have 
very small relative error of less than 0.2 per cent over the range of the scale factor 
1/1000 < a ^ 1 and the density parameter 0.2 < ^ 1 for a flat cosmology, and less 
than 0.4 per cent over the range 0.2 < ilm < 4 for a cosmology without a cosmological 
constant. 
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1 INTRODUCTION 



where, 



■ iHeathI j 19771 ) has shown that the growth function in a dust 
cosmology can be written as 
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H{a) = y^fln.a~'i + {1 - - !^A)a-2 + Qa, 



(1) 

(2) 
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where a is the scale factor of the universe, normalised to 
unity at present, Qm and Qa is the density parameters of 
dust matter and a cosmological constant, respectively. 

Although a compact expres sion using the incomp lete 
beta function has been shown bv lSildhauer et ali ()l992h . so 
far, no analytic solution of Di{a) has been presented for 
SIa 7^ 0. Here we restrict ourselves to the case ilm + i^A ~ 1, 
and present an approximate formula in a simple algebraic 

form. 

Following lEisensteinl (|l997l ). we adopt a normalization 
for Di{a) as 



da' 
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2 ' ' /o {a'H{a')) 
Then, in a fiat cosmology flm + = 1, our formula is 

(5) 



Note that our approximate formula is exact when fini ~ 1- 

2 APPROXIMATION 

We define 

F{a) = |n| r ^. (7) 

2 Jo {a'V^ma'-'i + 1 - Qm)) 

The power series expansion of F{a) around a = yields 

1 . 15 75 2 175 3 \ 

F(a) = a2 1 x+ -\ , (8) 

^ ' V 22 136 368 y ' ^ ' 

where x is defined in Eq. ©. After expanding F(a) up to 
0(x^^, we can obtain the Fade approximant to the following 
order: 

ll + M^+62^^+6^ 
1 + cia; + cix^ + ca^s 
where the numerical constants are determined as follows: 
619226202351 



D\(a) — aVT+x 



1 + 1.175j; + 0.3064a;2 + 0.005355a;^ 



1 + 1.857a; + 1.021^2 + 0.1530a;3 
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527102715964 ' 
12478731282519 
40730664415400 ' 

232758215919527 
43467765064114880' 
88964947071 
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Table 1. The maximal relative error AE (per cent) for the 
normalised growth function D{a) by our formula Eq. ^ for 

Cm + SIa = 1 



fin 



0.2 



0.3 



0.4 



0.5 



1 



maximal AE 0.20 0.03 < 0.01 < 0.01 



Table 3. The maximal relative error AE for the normalised 



growth function 


D(a) by 


our 


formula Eq. 


ig for 


rim 


0.2 


0.5 


1 2.0 


4.0 


maximal AE 


0.34 


0.02 


0.07 


0.36 



Table 2. A comparison of the relative error of Eq. II18I I and our 
Di(l) from Eq. igj. 



0.2 



0.3 



0.5 



0.9 



AE of Eq. (O 
AE of Eq. (0 



0.54 
0.19 



0.134 
< 0.01 



0.057 0.019 
< 0.01 < 0.01 



C2 = 



C3 = 



2445658735707 
2395921436200 ' 
17010766061223 
111170754639680 ' 
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We have checked the Taylor expansion and the Pade approx- 
imant using Maxima0 It is convenient to use a normalised 
growth function 

CW = §^, (10 

so that D(l) = l. 

In order to check the accuracy of our formula, we cal- 
culate the following relative error 



AE = 



jj-jappr _ jjnuu 
J^num 



X 100 (percent), 



(17) 



where D^pp"^ and D""™ represent the values of normalised 
growth functions calculated by using approximate formula 
and numerical method, respectively. 

Table [T] shows the maximal relative error AE for 
1/1000 ^ a ^ 1 in our method. It is apparent that our ap- 
proximate formula has sufficiently small uncertainties over 
the wide range of parameters 0.2 ^ flm ^ 1 and 1/1000 ^ 

The appr oximation of ICarroll et al.l (|l992l l. which was 
adopted from lLahav et aP ( 199ll ). is 
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(1 - Om) + (1 + 
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for flm + SIa = 1. 

A com parison of the relativ e error AE at a = 1 for both 
methods bv lCarroU et al. I (|1992| ) and us, is shown in Table [2] 
for f2m + f2A = 1. Our formula has generally smaller relative 
error in the range 0.2 < f2m < 1- 



in a slightly different context. Here, we show a simple and 
efficient approximate formula, directly as a function of a. In 
a dust cosmology with zero cosmological constant, SIa = 0, 
our formula is 



D 



i(a) = ay/l + y—^ 



1 + l. llSy + 0.1247y^ - 0.003893^ 
2.185y + 1.424J/2 + 0.24021/3 



where, 

y 



(19) 
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It is straightforward to derive Eq. H19|) from the Pade ap- 
proximant, in the same way as shown in the previous section. 
Table [3] shows the maximal relative error of the our for- 
mula for the normalised growth function D{a) in the range 
0.2 ^ f7m ^ 4 and 1/1000 S a ^ 1. 



4 CONCLUSION 

We have presented a simple approximation formula for the 
growth function in a spatially flat cosmology with dust and 
a cosmological constant. Our formula is written in terms of 
a rational function, and widely applicable over the range 
1/1000 < a ^ 1 and 0.2 < f2m ^ 1 with sufliciently smaU 
relative error of less than 0.2 per cent. 

We have also shown the approximate formula in a dust 
cosmology without a cosmological constant in terms of a 
rational function. The single rational function applies for 
all, open, closed and flat universes with relative error of less 
than 0.4 per cent over the range 1/1000 < a S 1 and 0.2 < 
fim < 4. 
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3 ZERO COSMOLOGICAL CONSTANT CASE 

In the Qa = case, the analytical form fo r Di(a) is widely 
studied (e.g., IWeinberg|[l97alHeathlll977l ). iBildhauer et al.l 
l|l992h also have shown the growth function as a function of 
the scale factor a in various Friedmann-Lemaitre universes. 
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